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ABSTRACT 

We  study  a discrete  model  based  on  the  observed  behavior  of  excitable 
media.  This  model  has  the  basic  properties  of  an  excitable  medium,  that  is, 
a threshold  phenomenon,  a refractory  period,  and  a globally  stable  rest  point. 
We  are  mainly  interested  in  two  dimensional  periodic  patterns.  We  characterize 
the  initial  conditions  which  lead  to  such  patterns,  by  introducing  a basic  in- 
variant, the  "winding  number  of  a continuous  cycle". 
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SIGNIFICANCE  AND  EXPLANATION 

Pattern  formation  in  living  organisms  is  a basic  problem  in  biology. 
Recently  chemists  have  been  studying  certain  (inorganic)  chemical  reactions 
which  lead  to  interesting  temporal  and  spatial  patterns,  in  an  effort  to 
understand  how  such  patterns  arise.  The  interaction  of  diffusion  and  chemical 
reaction  effects  seems  important.  In  this  paper  a discrete  mathematical  model 
is  analyzed  to  see  how  these  effects  can  lead  to  the  sorts  of  behavior  which 
have  been  seen  experimentally. 


The  responsibility  for.  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A COMBINATORIAL  PROBLEM  ARISING  IN  THE  STUDY 


OF  REACTION-DIFFUSION  EQUATIONS 
J.  Greenberg,  C.  Greene,  and  S.  Hastings 

The  problem  we  shall  consider  bears  a superficial  resemblance  to  the  well- 


known  "game"  of  Life,  as  devised  by  E.  J.  Conway  [1] , in  which  a set  of  simple 
rules  determines  the  step-by  step  evolution  of  certain  patterns  in  an  infinite 
planar  grid.  Our  problem  is  also  set  on  an  infinite  grid  of  square  "cells" 
and  proceeds  in  discrete  time  steps.  However  it  differs  from  its  predecessor 
in  having  a natural  physical  interpretation,  in  terms  of  reaction-diffusion 
processes.  These  are  of  current  interest  because  of  their  importance  in  a 
variety  of  biological  phenomena,  including  nerve  conduction  and  morphogenesis. 

A related  paper  [2]  continues  the  study  of  discrete  models  of  such  processes  which 
was  begun  in  [3] . However,  in  our  opinion,  the  combinatorial  aspects  of  the 
problem  have  sufficient  interest  to  warrant  separate  treatment. 

To  describe  our  process  we  label  cells  c = c^  ^ with  integer  coordinates 
(i/j)»  -0°  < i»j  < 00 , and  consider  an  infinite  sequence  t = 0,1, 2, 3,...  of 


discrete  time  steps.  To  each  triple  (i,j,t)  associate  an  integer  ufc  ., 

1 1 3 

called  the  "state"  of  cell  c.  . at  time  t.  These  integers  will  come  from 

i»3 

a fixed  finite  set  S = {0,1,2, ... ,n},  where  N ^ 2 . The  initial  states 

0 t 

u.  . are  chosen  arbitrarily  from  S.  Subsequent  states  u.  .,  t > 0,  are 
1* 3 i,3 

then  determined  inductively,  according  to  rules  A and  B below. 

The  inductive  procedure  to  be  described  requires  choosing,  initially,  a 


fixed  integer  K with 


1 < K < N/2  . 
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The  motivation  for  these  rules  is,  roughly,  that  excitation  "diffuses" 
from  an  excited  region  into  an  adjacent  resting  region,  but  not  into  a 
refractory  region.  Also,  once  a cell  is  excited,  its  state  evolves  according 
to  fixed  dynamics  with  no  diffusion  effects  from  neighboring  cells,  until 
it  returns  to  rest.  We  remark  that  a number  of  previous  authors,  starting 
with  Wiener  and  Rosenblueth  [4] , have  studied  similar  processes,  usually  on 
a computer  and  again  in  a biological  context.  In  [2]  it  is  shown  how  these 
rules  are  related  to  a certain  singular  limit  of  some  widely  studied  continuous 
models  of  reaction-diffusion  processes. 

The  problem,  broadly,  is  to  describe  how  a given  initial  pattern 


P = ((i,j,u°  .),  < i,j  < “>}  evolves  as  t increases.  In  particular, 

u i , j 


what  sorts  of  patterns  can  develop,  and  will  the  process  continue  indefinitely 
without  all  cells  returning  eventually  and  permanently  to  rest.  In  [3]  it  is 
observed  that  for  N = 2 and  K = 1 there  is  a complete  solution,  provided 
only  that  the  number  of  non-zero  cells  at  t = 0 is  finite.  There  are  two 
possibilities. 


I.  The  pattern  dies  out.  By  this  we  mean  that  for  any  (i,j)  there  is  a 


T.  , such  that  u.  . = 0 if  t > T.  ..  Equivalently, 
i*3  - 1*3 
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(2) 


d(m,n)  = min{ |ra  - n|,  N + 1 - |m  - n | } 


for  any  m,n  in  S.  Equivalently,  identify  each  k in  S with  the  point 


k = e 

on  the  unit  circle  C in  the  complex  plane.  Then 

d(m,n)  = N (shorter  distance  from  m to  n on  0}  . 

2v 

Observe  that  for  any  cell  c.  . and  any  t > 0, 

1,3  - 

(3)  d(ufc  .,ut+l)  < 1 . 

1,3  1,3  - 

Assume  that  N 3 and  let 

(4)  L = min(K  + 1,  N-  / 1 } . 

4 

t t 

Theorem  1.  If  there  is  a t_  > 0 such  that  d (u  . ,u..  ..)  < L whenever 

0 - u , 3 1,3  

cij  and  cit  j.  are  adjacent,  then  the  pattern  dies  out.  (Recall  that 

c.  -i  and  ci  . , are  adjacent  if  |i  - i'|  + |j  - j'l  = 1.) 

In  particular,  if  the  process  is  persistent,  then  there  must  be  adjacent 

cells  c and  c.,  such  that  d (u°  .,u°,  > L.  In  fact,  this  can 

i,3  1,3  i»3i»3— 

be  strengthened  a bit. 

Theorem  2.  If  the  process  persists,  then  there  is  a fixed  pair  of  adjacent 


cells  c.  . and  c..  ..  such  that 

i,3  1,3  


iL 


for  all  t > 0. 


N + 1 

If  L < — - — , then  eventual] y an  even  stronger  discontinuity  must 
develop. 

Theorem  3 . If  the  process  persists,  then  for  any  sufficiently  large  t there 

is  a pair  of  adjacent  cells  c.  , and  c,.  ..  such  that  d (ufc  .iifc  ^ > N + 1 

i,3  i .3  i»  j i • f j • ' - — 4 — 
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An  estimate  will  be  given  for  when  (at  what  t)  this  inequality  must 


hold . 

In  order  to  give  sufficient  conditions  for  persistence  we  introduce  the 
concept  of  a cycle. 

12  3 M m+1 

Def.  A cycle  is  an  ordered  (M  + 1) -tuple  C = (c  ,c  ,c  , ,c  ,c  ) of 

1 M m+1  1 

cells  such  that  c , ...,c  are  distinct,  c = c , and  c is  adjacent 

to  c1+1  for  i < i < M. 


(FIG.  3) 

Def.  A cycle  C is  said  to  be  continuous  at  time  t if 

d(ut,u!:  ) < K for  1 < i < M , 
l l+l  — — — 

where  u^  is  the  state  of  cell  c1  at  time  t. 

For  such  a cycle  we  then  define  a "winding  number"  at  time  t.  For  this 

purpose  recall  the  previous  identification  of  the  states  k = 0,1,2, ,N 


with  the  points  k = e 


on  the  unit  circle.  If  m,n  e S,  let  mn 


denote  the  shorter  directed  arc  from  m to  n.  If  both  arcs  from  m to  n 
are  of  the  same  length,  let  mn  be  the  arc  connecting  m to  n in  the 
counter-clockwise  direction.  Then,  for  an  ordered  pair  (m,n)  of  integers 


in  S,  let 


o(m,n)  = 


d(m,n)  if  mn  connects  m to  n in  the  counter-clockwise  direction 
-d(m,n)  otherwise  . 


The  winding  number  of  a continuous  cycle  C at  time  t is  then  defined  by 

Vc>  * FTT  • 

1*1 

It  is  not  hard  to  show  that  (Cl  is  an  integer,  and  represents  the 

net  number  of  times  the  unit  circle  is  traversed  in  the  counter-clockwise  direc- 


tion by  the  points  u^^  as  i runs  from  1 to  M + 1.  We 
and  sufficient  condition  for  persistence. 


now  give  a necessary 
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a continuous  cycle  C at  time  T such  that  WT (O  * 0. 

Obviously  this  includes  a sufficient  condition  for  persistence  which 
can  be  checked  at  t = 0. 

It  is  desirable  to  find  an  upper  bound  for  the  smallest  T satisfying 
the  conditions  in  Theorem  4.  I-  can  be  shown  by  example  that  T may  be 
arbitrarily  large  if  the  size  of  the  initial  non-zero  set  is  not  restricted. 
Our  result  in  this  direction  is  probably  not  the  best  possible. 

Theorem  5.  Let  R^  be  a "diamond"  shaped  set  of  the  form 

|i|  + | j | < m + 1} 


R.  = {c.  . 

1 1 » J 


for  some  m,  and  suppose  that 

u°.=0  if  I i I + I j I > m . 

1,3  1 1 11 

Then  the  pattern  persists  if  and  only  if  there  is  a continuous  cycle  with  non- 
zero winding  number  at  time  T = N(R^)  • K,  where  n(R^)  is  the  number  of 
cells  in  R^. 

Our  final  result  is  proved  in  almost  the  same  way  as  Theorem  5. 

* 

Theorem  6.  The  process  is  eventually  periodic  in  any  bounded  region. 

In  other  words,  there  is  a p such  that  for  each  m there  is  a T with 

in 

t+p  t 

u . . = u . . 

lO  ir] 

i I + I j|  1 ra. 

One  can  extend  these  results  in  several  directions.  For  instance,  arrays 

of  cells  in  more  than  two  dimensions,  or  with  non-rectangular  geometry  could 

be  considered.  Other  definitions  of  "adjacent",  or  "neighboring"  cells  might 

be  used.  For  example  in  the  plane  with  square  cells  we  could  say  that  c.  . 

1,3 

and  c.,  are  neighbors  if  l<|i-i'|+|j-j'|i2*ln  any  of  these  cases 

i • 3 


whenever  t > T and 
— M 
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Theorem  4 goes  over  without  change.  However  Theorem  1 may  need  modification. 
The  alternative  definition  of  adjacent  cells  given  above  requires  the  number 


N + 1 


N + 1 

instead  of  — : — in  Theorem  1.  On  the  other  hand,  with  rectangular 


3 4 

geometry  and  the  definition  of  two  cells  as  adjacent  if  they  have  common  faces, 

Theorems  1-6  are  essentially  unchanged  in  higher  dimensions. 

Proofs : A basic  observation  is  that  discontinuities  in  a cycle  do  not  appear 
spontaneously.  This  is  implied  by  the  following  result. 

Lemma  1 . Suppose  that  c and  d are  adjacent  cells  with  states  u*"  and  vfc 
at  time  t.  If 


(5) 


for  some  tQ  >_  0,  then 


d (u  °,v  °)  < K 


t t 

d(ut,vt)  < max{d(u  °,v  °)  , 1} 


for  all  t >_  t . 

Corollary . If  a cycle  is  continuous  at  t^,  then  it  is  continuous  for  all 

i V 

t t 

Proof  of  Lemma  1 . If  u ^ 0 and  v ^ 0,  then  rule  B,  together  with 
(2) , implies  that 


(6) 


t +1  t +1  t 

. , 0 0 . . , 0 C> 

d (u  ,v  ) = d (u  , v ) . 


We  can  therefore  assume  that  at  least  one  cell,  say  c,  is  in  the  resting 

t t 

state  at  t = t , i.e.  that  u =0.  If,  in  addition,  v =0  then  rule  B 

V1  V1 

implies  that  d(u  ,v  ) < 1. 


Next  suppose  that  u 


0,  1 v K.  Again  use  rule  B,  to  conclude 

t. 


V1 

that  u =1  and  again  (6)  holds.  Finally,  if  K+l<^v°<^N  then 
(1),  (2),  and  (5)  imply  that  v^>N-K+l>  N/2.  But  from  this  it  follows 
that 
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t +1  t t t 

o 0.  0 0.  . 

d (v  ,u  ) = d (v  ,u  ) - 1 . 

tQ+i  to+1  fc0  fco 

From  (3)  and  the  triangle  inequality  we  get  d(u  ,v  ) <^  d (u  ,v  ), 


completing  the  proof  of  Lemma  1. 

It  turns  out  that  Theorem  4 is  the  key  result  so  we  prove  it  first. 

Proof  of  Theorem  4.  We  begin  by  showing  that  if  C is  a continuous  cycle  at 
t * and  hence  for  t t , then  W (C)  - W (C)  for  t >_  t . It  suffices 

w U U Lq  U 

to  show  that  W (C)  = W (C) . 

fc0+1  fc0 

For  each  pair  of  integers  j and  k with  1 < j < k M , let 


2t(j'k)  = l.  °{Vui+i) 

1=3 


In  particular,  W (C)  = ■ ■ ■ Q (1,M  + 1). 

t N + 1 t 

Lemma  2.  Suppose  that  j and  k are  integers  with  l£j<k£M+l  and  that 

(7)  u*"  ^ 0,  u^  * 0,  and  ut  = 0 if  j<i<k. 

3 k l 


Then  Qt(j,k)  =Qt+1(3'k)* 

Proof.  If  k = j + 1,  then 

Qt(j,k)  = 0 (Uj 'u^)  = 0 (Uj+1 >u^+1)  = Qt+1 ( j ,k) 

/\ 

s s 

because  both  u . and  u move  one  step  counter-clockwise  around  the  unit 
3 k 

circle  as  s goes  from  t to  t + 1. 

Now  suppose  that  k > j + 2.  From  (7)  we  see  that 


0 or  1 if  j<i<k. 


„ , t+1  t+1,  t+1  t+1  . . . 

Hence  o(u.  ,u.  „)=u.  ,-u.  if  n < l < k - 1 and  so 
l l+l  l+l  l 


Qt+l(j'k)  = + ._I+1  Kti  - ur1)  - a(u£,u*+1>  . 


where  the  summation  term  on  the  right  is  not  present  if  k = j + 2.  If 
k > j + 2 then  the  summation  term  collapses.  We  conclude  that  for  any 
k 1 j + 2, 


<8i  Afi'j.w  ■ + vi  - vi  * “vX*1'  • 

From  (7)  it  follows  that  ufc  ^ and  u.t+^  lie  in  the  set 

J k 

[N  - K + 2,N]  U [0,K  + 1]  . 

t+1  . t+1 

Since  and  u ^ are  0 or  it  ls  easily  seen  that 

t+1  t+1.  t+1  , t+1 

° uj  'uj+l}  - uj+l  = 0(uj  '0) 


„ . t+1  t+1,  j t+1  t+l 

° uk-l'uk  ) + Uk-1  = °(0'uk  > 


and  from  (8) , 

Qt+l(j,k)  = °<uj+1*°>  + °(0,u^+1) 

= o(u^,0)  + a(0,uS 
0 k 

= Qt(j»k)  , 


where  we  again  use  (7).  This  proves  Lenina  2. 

The  first  part  of  Theorem  4 then  follows  quickly  by  letting 
i-2  < ■*’2  < * * * < ’’•p  t*1086  i such  that  u^  ^ 0,  and  observing  that 


VC)  * j =t‘V  W • 


where  we  set  i ^ = T^e  desired  result  follows  by  applying  Lemma  2 to 

j = i , k = i for  1 < m < p. 

Remark.  An  alternative  proof  that  the  winding  number  of  a continuous  cycle 
is  constant  proceeds  as  in  the  following  sketch:  The  states  of  the  cells 
of  C at  time  t can  be  used  to  define  a continuous  map  F : C -+  C of  the 
unit  circle  into  itself.  To  do  this,  identify  the  cells  c ,...,c  of  C 


/s. 

with  the  points  c3  = e 


2iri  /N 

on  C.  Map  c3  into  u_..  Extend  this  map 


to  a continuous  one  from  C to  C by  interpolation,  using  shortest  arcs 
along  C and  taking  the  counter-clockwise  arc  in  case  of  ties.  This  defines 
F^,  and  similarly  one  can  define  , using  the  states  at  time  t + 1. 

These  two  maps  are  easily  seen  to  be  homotopic,  and  hence  have  the  same 
winding  number. 

We  next  consider  the  second  part  of  Theorem  4,  namely,  that  if  a pattern 

persists,  then  eventually  there  must  be  a continuous  cycle  with  non-zero 

winding  number.  To  prove  this  some  additional  concepts  are  helpful. 

1 M i 

Definition.  A path  is  an  M-tuple  P ® c ,...,c  of  cells  such  that  c is 


adjacent  to  c 


i+1 


for  1 < i < M - 1. 


Def.  A path  P is  said  to  be  continuous  at  time  t if  d(ut,ut.,)  < K 
i l+l  — 

for  1 < i < M - 1,  where  ufc  is  the  state  of  c1  at  time  t. 

— — i 

By  Lemma  1,  if  P is  continuous  at  tQ,  then  it  is  continuous  for  all 

We  shall  say  that  a cell  c is  external  to  the  pattern  at  time  t if 

there  is  a rectangle  R in  the  plane  such  that  c is  outside  R but  all 

cells  which  have  non-zero  states  at  time  t lie  inside  R. 

Under  the  assumption  that  there  are  no  continuous  cycles  with  non-zero 

winding  number  for  any  t > 0,  we  can  define  the  potential  of  a cell  c.  . 

i,l 

at  time  t for  any  cell  which  is  connected  to  an  external  cell  by  some 


1 M M 

continuous  path  at  time  t.  Let  the  path  be  c ,...,c  , where  c 
and  c3-  is  external.  Then  we  set 


c . . 
i,l 


M-i 

hfc(i,j)  = l °<VUi+l) 


k=l 


This  will  be  the  same  for  any  continuous  path  connecting  c„  to  any  external 
cell,  for  otherwise  we  could  find  a continuous  cycle  with  non-zero  winding 
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number.  All  cells  external  to  the  initial  pattern  have  a potential  for  all  t. 

Also,  if  h (i,j)  is  defined,  then  so  is  h (i,j)  for  t >_  t , and 
0 to 

h (t,j)  >_  h (i,j),  t t . Among  those  cells  which  are  not  external  to  the 
0 

initial  pattern  (a  finite  number),  we  can  allow  the  possibility  that  some  may 

never  have  a potential  defined.  In  any  case  there  must  be  a t^  >_  0 such 

that  no  cells  have  a potential  defined  for  the  first  time  at  some  t >_  t^. 

Let  c.  . be  any  cell  with  a potential  defined  for  t > t . Suppose 
l / D 0 

that  for  some  t >_  t^,  u^  _.  = 0 and  1 <_  u^,  K for  an  adjacent  cell 

c ■ i .it-  Then  c.,  has  a defined  potential  at  time  t which  must  be 

1 » j i » j 


by 


higher  than  that  of  c.  .,  since  we  can  connect  c..  ..  to  the  outside 

1*3 

a path  going  through  c.  .. 

1 * j 

Let  A'  be  the  highest  potential  of  any  cell  at  time  t^.  Let  A be 

the  next  higher  multiple  of  N + 1.  We  claim  that  no  cell  can  achieve  a 

potential  greater  than  A.  If  not,  let  c be  the  first  cell  to  achieve  a 

k,Z  t -1 

potential  of  A + 1,  and  suppose  this  occurs  at  t = t . Then  u =0, 

i K f jC 

and  c has  a neighbor  c , , which  is  excited  at  time  t - 1,  so  that 

K $ X,  K t a,  ±. 

"r1 

1 £ u , , K.  But  then  c , must  have  a potential  greater  than  A at 

K t jL  K t jC 

t = t^  - 1,  which  is  a contradiction. 

But  this  proves  that  among  all  those  cells  with  potential  A'  at  time 
tQ,  none  can  ever  become  excited  again,  contradicting  the  persistence  of 
the  pattern. 

Proofs  of  Theorems  1,  2,  and  3.  Theorems  1,  2,  and  3 all  follow  from  Theorem  4. 

Notice  that  as  a consequence  of  Lemma  1,  Theorems  1 and  2 are  corollaries  of 

N + 1 

Theorem  3.  To  prove  Theorem  3,  assume  that  i = K + 1 < — ^ — and  choose 

tQ  large  enough  to  insure  that  at  tQ  there  is  a continuous  cycle 
1 M 

C = (c  ,...,c  ) with  W ( C ) * 0.  This  is  possible  by  Theorem  4.  We  shall 
show  that 
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t t 

, , 0 0 , . N + 1 

d (u.  . .u.  . , ) > — - — 

1 - 4 

for  some  pair  of  adjacent  cells  c.  . and  c., 

1 » D i » D 

For  this  purpose  we  use  a different  kind  of  continuity  for  a cycle.  We 
say  that  a cycle  E = (e\...,e^)  is  mildly  continuous  at  t^  if 


d(u.  ,u  ) < 
l l+l  — 


N + 1 


for  1 < i < Q,  where  u^  is  the  state  of  e at  t . (We  shall  only  be 
concerned  with  states  at  t = t^,  and  so  we  suppress  the  time  in  this 
notation. ) 

For  any  mildly  continuous  cycle  E,  the  winding  number  W (E)  is  defined 

t0 

just  as  before.  We  are  assuming  that  W (E)  ^ 0. 

t0 

t t 

Let  r = sup{d(u.  . ,u.,  ,,)/c.  . is  adjacent  to  c.,  .,}  and  suppose 
1 * 3 1 » 3 1 » 3 i i] 


that  r < 


N + 1 


Lemma  4.  There  is  a mildly  continuous  cycle  E at  tfi  which  consists  of 
exactly  four  cells  and  has  non-zero  winding  number. 

Proof.  Consider  C as  an  ordered  M-tuple  of  closed  squares  in  the  plane. 

Let  y be  the  Jordan  curve  obtained  by  joining  the  centers  of  consecutive 
squares  of  C.  Suppose  that  the  inside  of  y contains  squares  which  are  not 
in  C.  Such  cells  we  call  interior  to  C.  Then  there  must  be  three  consecu- 
tive cells  c1,c1+^’,c1+^  in  C which  form  three  quarters  of  a square  of 
four  cells  such  that  the  fourth  cell  c*  in  this  square  lies  inside  y. 

(FIG.  4) 

Since  r < N-  ~ - , the  cycle  obtained  by  replacing  c1+*  in  C by  c*  is 

still  mildly  continuous.  Furthermore, 

0(ui'ui+2)  " °(ui'ui+l)  + a(ui+l'Ui+2) 

= 0(ui<u*)  + 0(u*,Ui+2) 
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where  u*  is  the  state  of  c*  at  tQ.  Therefore  the  new  cycle  C*  has  the 
same  (non-zero)  winding  number  as  C.  However,  c*  has  one  less  interior 
cell  than  C. 

Continuing  this  shrinking  process,  we  see  that  there  must  be  a mildly 

continuous  cycle  C 3 at  tQ  with  non-zero  winding  number  and  no  interior 

cells.  Again  there  must  be  consecutive  cells  c1,c1+*,c1+2  of  C * which 

comprise  three  quarters  of  a square  of  cells,  and  now  these  can  be  chosen  so 

* 1 

that  the  fourth  cell  c in  this  square  is  also  a cell  of  C . We  can  renumber 

A “I 

the  cycle  so  that  c = c for  some  j > i + 2. 

If  j = i + 3,  then  there  are  two  cases.  The  cycle  C = c1,c1+^’,c1+2,c*+3 
may  have  non-zero  winding  number,  in  which  case  we  are  done.  If,  on  the  other 

a i+i  i+2  1 

hand,  C has  winding  number  0,  then  omit  c and  c in  C and 

there  results  a cycle  C2  with  W ( C 2)  = w (C1) . Thus  C1  has  been  reduced 

0 fcc 

to  a still  smaller  cycle. 

Next  suppose  that  j > i + 3.  Proceed  according  to  whether  the  cycle 

D = (c1  2,c1+3 CJ)  has  non-zero  winding  number  or  not.  If  W ( D ) ^ 0, 

replace  C with  C.  If  w ( D ) = 0,  eliminate  c1+2,...,c^  from  . 

t0 

It  is  thus  clear  that  we  arrive  eventually  at  a cycle  with  the  desired 
properties  at  t , proving  Lemma  4. 

But  from  this  result  we  obtain  a contradiction.  Since 


d(u. ,u.  ) < r < 

l l+l  — 


N + 1 


and  assuming  as  we  may  that  u^  = 0,  we  have  the  inequalities 

1 < u_  < r,  u,  < 2r,  u.  < 3r  . 

— 2 — 3 — 4 — 

On  the  other  hand,  u^  >_  N - r + 1.  This  gives 

N + 1 

3r  > N - r + 1,  or  r>  — ; — , 

— — 4 

which  contradiction  proves  Theorems  1,  2,  and  3. 
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Proof  of  Theorems  5 and  6.  Let  be  the  diamond  shaped  region 


R = (c.  . 

0 1,3 


I i I + I j I 1 m)  * 


where  m is  chosen  so  large  that  RQ  contains  the  entire  non-zero  set  at 
t = 0.  Also,  let 


R = (c . 
n 1,3 


| x | + |j|  _<  m + n) 


for  n = 1,2,3,...  . 

Lemma  5.  Within  any  given  R^,  n ^ 1,  the  process  proceeds  independently 

from  the  outside  of  R . More  precisely,  if  c . . e Rn  and  u.*3.  = 1,  then 

n — 1,3  1,3  

2 < u.?  < K + 1 for  some  cell  c.  , adjacent  to  c.  ..  and  also  contained 

1 , j 1 » j 1 » 3 

in  R . 

— n 

This  result  implies  that  if  j is  a second  process,  obeying  the  rules 

A and  B if  c . . e R°  but  following  the  rule  u!”  . = 0 if  c . . It  R°, 

t*3  i»D  i, j 

then  ufc  . = ufc  . in  R°.  The  process  u*”  . is  clearly  eventually  periodic, 
i,3  i*3  i*3 

which  leads  to  our  theorem. 

Proof  of  Lemma  5.  Suppose,  then,  that  the  Lemma  is  false,  and  let  t^  be  the 

first  time  where  a cell  in  some  R , n > 1,  is  excited  by  a cell  outside  R , 

n — n 

C 

and  not  simultaneously  excited  by  a cell  in  R . Thus  there  is  a cell  c.  . = c 

t t t n 1,3 

in  R with  u.°.  = 11  0 = 1,  while  2 < u,  < K + 1 for  some  adjacent  cell 
n i,  j 0 — 1 — 

, n+1 

c in  R . Furthermore, 


(8) 


u2u  * [ 2, K + 1] 


2 n . 0 

for  any  cell  c in  R which  is  ad3acent  to  c . 

However  u^  - u^  = 0,  so  Lemma  1 implies  that  u^°  = 2.  Also,  c3  must 
in  turn  have  been  excited  by  a cell  c3  6 Rn.  This  follows  from  our  hypotheses 
on  tQ,  since  cells  of  Rn+1  which  are  not  in  Rn  cannot  be  adjacent  to  each 
other.  Then  u^0  - 3.  Also,  c°,  c1  and  c3  form  three  cells  out  of  a 
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2 

square  of  four  cells.  Let  c be  the  fourth  cell  in  this  square.  Then 

c2  e Rn  and  c2  is  adjacent  to  c°  and  c3.  Since  u ° = 3,  urt°  = 1,  and 
t t t t t 3 

d (u0  ,u2  ) 1 d^u3  ,U2  ) i it  is  seen  that  u ° = 2, 


3 '"2 

of  (8) . This  proves  the  Lemma. 


3 ' 0 

a contradiction 


Theorem  6 is  immediate  since  the  process  u!"  . , which  is  on  a finite 

1 * J 

grid,  is  clearly  periodic.  Theorem  5 also  follows  easily,  since  u*”  . can 

1*3 

only  be  persistent  if  there  is  a cycle  c1,...,cN+1  of  cells  in  R such 

M+l 

that  1 Uq  K,  and  each  c^  excites  cell  c^  3 within  a time  K from 

when  c1  was  excited  by  c1  1,  for  2 < i < N.  The  existence  of  such  a 


cycle  follows  easily  for  a finite  grid. 
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